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^ Abstract 

In this work we perform a first-principles study of the adsorption properties of an In adatom 
deposited on 1.75 monolayers (ML) InAs, forming a wetting layer on GaAs(OOl) with the 02(2 x 4) 
or /?2(2x4) reconstruction. The structural properties of these reconstructions have been studied: we 
determine the equilibrium geometry of the surfaces and their stability for various growth conditions. 

s 

^ We have then carried out a detailed study of the potential energy surface (PES) for an In adsorbate, 

^ finding the minima and the saddle points. The main characteristics of the PES and the bonding 

configurations of the In adatom on the surface are analyzed by comparing with analogous studies 
^ reported in the literature, trying to extract the effects due to: (i) the compressive strain to which 

the InAs adlayer is subjected, (ii) the particular surface reconstruction, and (iii) the wetting layer 
composition. We found that, in general, stable adsorption sites are located at: (i) locations besides 
the As in-dimers, (ii) positions bridging two As in-dimers, (iii) between two adjacent ad-dimers 
'""I (only in (52), and (iv) locations bridging two As ad-dimers. We find also other shallower adsorption 

sites which are more reconstruction specific due to the lower symmetry of the a2 reconstruction 

00 

compared to the /?2 reconstruction. 
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I. INTRODUCTION 



Thanks to the strong charge carrier confinement, QDs have atomic-hke properties that 
could be useful for applications in optical and optoelectronic devices, quantum computing, 
and information storage. One of the most challenging problems for quantum dot formation, 
is the control of their shape, composition and density. From this point of view the abundant 
experimental knowledge has not yet produced a sufficiently deep understanding of the physics 
of formation of quantum dots that could guide the way to control the growth process. In 
this situation first-principles simulations can be helpful in sheding light on the atomistic 
mechanisms that lead to dot nucleation at surfaces. 

In the present paper, we address one of the most-studied systems: InAs quantum dots 
grown on a GaAs(OOl) substrate, where the lattice mismatch between InAs and GaAs is 
of the order of 7%. Strain relaxation at the surface acts as a driving force for a self- 
assembly process, the so called Stranski-Krastanov growth mode [1], with the deposition 
of an atomically thin InGaAs wetting layer (WL) and subsequent surface mass transport 
accompanied by a 2D to 3D growth transition. 

We focus our attention on the structural properties of the WL on which quantum dots 
nucleate, and on the interaction between a single In adatom and the surface, in order to 
find where the diffusion barriers are located and what the possible adsorption sites on the 
surface function of the wetting layer reconstruction. 

We have chosen to investigate mainly the 0:2(2 x 4) reconstruction of the (001) surface 
because of its stability under In-rich conditions and under compressive strain |2l |3[ E] which 
is the case for the InAs WL grown on GaAs. Moreover, a (2 x 4) reconstruction was found 
[3 El [71 E] to occur on the (001) wetting layer at the onset of dot nucleation in a In-rich WL 
condition. We have modeled a WL at a InAs coverage 9 = 1.75 monolayer (ML) which is 
close to the critical value corresponding to the 2D to 3D transition [9J. The structural and 
geometrical properties of this surface have been investigated. Using the 02 reconstruction 
we have calculated the PES for a single In adatom, with the aim to understand what are the 
implications of the WL morphology (strain and composition) on the adsorption and surface 
mobility of the deposited In. To this end a comprehensive comparison with results reported 
previously in the literature for Ga and In adatoms and other surface reconstructions is given. 
In this context, a detailed analysis of the differences and analogies between the 02 and the 
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reconstructions is pursued. 

In the literature similar calculations have been reported for GaAs [100] homoepitaxy [TU] . 
The issue of the WL composition for In adsorption has been addressed only by Penev et. 
al. ^lll |12], but for a much lower In coverage of 6' = 0.66 ML. 

The paper is organized as follows: in section II we detail our theoretical approach. Section 
III reports the obtained results for: (i) the WL surface reconstruction and energy; (ii) the 
PES of a single In adatom on the 0:2 and the (32 reconstructed WL; (iii) the determination 
of the barrier positions and heights between each couple of adsorption sites; (iv) a detailed 
study of the atomic bonding of the In adatom to the WL. In section IV our results are 
compared with those regarding other surface reconstructions in order to extract a rationale 
for the dependence on In adsorption and diffusion on the local structure of the (001) surface. 
Finally, a summary is given at the end of the paper. 

II. THEORETICAL METHOD 

We performed first principles calculations within Density Functional Theory in the local 
density approximation (DFT-LDA) [13] using the exchange and correlation potential of 
Ceperley and Alder [13] as parametrized by Perdew and Zunger [T^. For this kind of 
calculations both LDA [21 [ISl [13 CSl CS] and GGA [H [HI EOl ED [22] have been used in the 
literature. It has been demonstrated [231 [23] that both methods give qualitatively the same 
picture for the description of adatom adsorption on a surface, although LDA usually gives 
larger values for binding energies than GGA [25j. Only the details can be different: LDA 
tends to overestimate the binding energy but for lattice constants, elastic moduli and surface 
energies the LDA calculated values are generally closer to the experimentally determined 
ones. We have used norm-conserving pseudopotentials treating the outermost s- and p-shells 
of Ga, In and As as valence electrons, and the electronic wave functions were expanded in 
plane-waves, with a 18 Ry kinetic energy cutoff. The energy cutoff has been tested in order 
to reach convergence for the lattice bulk properties of GaAs and InAs. The core corrected 
atomic pseudopotentials have been tested on bulk Ga, In and As, where the determined 
equilibrium configurations and elastic moduli compare well with the experimental data. 

The equilibrium lattice parameters obtained for the GaAs and InAs bulk phases are 
ao = 5.609 A and oq = 5.906 A, which are slightly smaller than the experimental ones. 
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Thus our calculated lattice mismatch is 5.4%, smaller than the experimental value. This is 
going to underestimate the WL strain. In sec. IV we will compare our results with results 
[TUt \T2[ [TB] where this effect was instead overestimated. 

Starting from the GaAs structure, we have set up the (001) oriented supercell containing 
4 layers of GaAs, covered with 1.75 layers of InAs arranged according to the 02(2 x 4) surface 
reconstruction. The lower layer of Ga atoms is kept fixed during the cell relaxation, in order 
to mimic the constraint due to the underlying semi-infinite bulk, and it is passivated with 
pseudo hydrogen atoms of 1.25 electron charge. The slab is repeated along the (001) direction 
with a periodicity 5ao and a separation of about 10 A of vacuum. As for the number of 
GaAs layers, we have found that only two layers suffice to obtain a correct description of the 
PES, minima and maxima positions. The refinement using two more GaAs layers produces 
a difference in the barrier heights all within 40 meV. Brillouin-zone (BZ) integration was 
carried out using a set of special fc-points equivalent to 16 points in the 1x1 surface BZ. A 
smearing of 0.02 Ry has been used in order to account for the possible metalization of the 
surface electronic structure. 

In the ^2(2 X 4) reconstruction the last complete atomic layer is constituted by As atoms. 
A 6 = 0.75 cation layer is deposited over it. This layer is then terminated with one As dimer 
(As ad-dimer) on top (see Fig. [T]). The uncovered As on the last complete layer dimerizes 
along the [110] direction (As in-dimer). The /32(2 x 4) reconstruction is similar but presents 
one additional As ad-dimer bridging the two remaining In rows on the top of the surface (see 
Fig. [1]). It has a larger symmetry than the 02 reconstruction, being mirror symmetric with 
respect to the (110) plane passing through the As in-dimers. All the examined structures 
have been relaxed in order to find the equilibrium geometries, until all the forces acting on 
the atoms were less than 2.5 ■ 10"'^ eV/A. 

The PES of an In adatom is calculated by relaxing the adatom z component together 
with all the surface degrees of freedom while the in-plane x, y coordinates are kept fixed. 
We have set up a grid of 8 x 4 points, corresponding to a step of 2 A, along the (110) 
and (110) directions respectively, and for each point we have found the minimum energy 
configuration for the adsorbate plus surface. Then we have interpolated the grid data with a 
bi-cubic spline algorithm in order to find the positions of the minima and the saddle points. 
The exact values of the minima have then been determined by further relaxing all the three 
coordinates of the adatom together with the surface. To minimize the adsorbate interaction 
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FIG. 1: (Color online) Representation of the 02(2 x 4) (left) and (32{2 x 4) (right) surface recon- 
structions. The shown geometries are those optimized. It is evident the bending of the the In-In 
bond in the top layer in the 02 reconstruction. Dark balls and light balls represent In and As ions 
respectively. 

we have doubled the surface unit cell along the (110) direction. 

Since one of the objectives of this work is ultimately to study In diffusion on the WL 
within the framework of the Transition State Theory [26 1 127 1 128]. the minimum energy paths 
(MEP) and the height of each energy barrier between every couple of adsorption sites have 
to be calculated, leading to a more detailed knowledge of the PES. This task has been 
accomplished by using the Nudged Elastic Band method (NEB). This method is able to 
find the MEP between two adsorption sites, by simulating a string of replicas of the system, 
where the different images are one linked to the other by springs |29t I5U] . By minimizing the 
energy associated to the path, an accurate description of the MEP and thus of the saddle 
point is obtained. 

All the ab-initio calculations were performed using the ESPRESSO simulation package 
ETl. 



III. RESULTS 



A. Structure of (001) (2 x 4) surface reconstructions 

We have calculated the optimized geometries and the formation energies of the 02(2 x 4) 
and /?2(2 x 4) reconstructions of the bare GaAs surface and of GaAs covered with a ^ = 1.75 
ML InAs layer. The surface energies vs. the As chemical potential ^as (whose value is 
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related to the growth conditions such as growth temperature, cation/anion flux ratios) are 
calculated through the expression [32] 

7f = ^ [Etot - ^'Ga/UGaAs - ^'inAHnAs) - /"As , (1) 

(/^inAs = 0, A^in = for purc GaAs). The obtained surface energies are reported in Fig. |2} 
Here i^tot is the total energy of the slab, yUxx is the chemical potential of the bulk phase of 
material xx, A^y is the number of atoms of element y in the supercell and A is the surface unit 
cell area. A set of calculations is performed for a slab with two hydrogen passivated surfaces, 
in order to subtract the contribution of the bottom surfaces of the reconstructed slabs to 
be investigated, hence i^tot is the net contribution of the reconstructed slab. Elemental As, 
Ga and In are calculated using their ground state structures: respectively the rhombohedral 
A7 structure for As, a-Ga and bet In. All these calculations have been done using high 
convergence standards. We see that the range of stability of the 02 phase, when GaAs is 
covered with the InAs ad-layer, is larger than that of the bare GaAs surface. This means 
that the formation of the 02 phase will be favored over the (32 phase over a larger range 
of As flux values (or surface growth conditions), when the InAs WL is formed. This is in 
agreement with previous calculations showing a stabilization of the 02 phase for InAs under 
isotropic compressive strain [21 |TT] which is the condition of the InAs layer grown on GaAs. 

Regarding the atomic structures, the comparison of the distances between the atomic 
planes, in the InAs covered GaAs case, shows that the values are systematically larger than 
those obtained for the pure GaAs surface by about 10 %, owing to the stress originated by 
the InAs/GaAs lattice mismatch. In particular, we remark the buckling of the bond between 
the In atoms in the uncompleted In layer found for the 02 reconstruction (see Fig. Q due to 
the missing (with respect to the [32 reconstruction) second As ad-dimer. It was shown ([2]) 
that this buckling plays a role in the stabilization of the 02 phase. 

B. The potential energy surface of the single In adatom 

In the following we describe our calculations for both the 02 and /?2 surface reconstruc- 
tions, with an accurate description of the adsorption sites and hopping barriers. 
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FIG. 2: (Color online) Surface energy phase diagram for the 02(2 x 4) and /?2(2 x 4) reconstructions 
of InAs/GaAs(001) and GaAs(OOl). 

1. a2(2 X 4) surface reconstruction 

For this surface reconstruction we have found eleven adsorption sites for the In adatom, 
and the related saddle points. The energies of the minima and of the saddle points are given 
in Tables [T] and [11} respectively, and a graphical representation of the PES is reported in 
Fig. [3j In the Tables, the energies are relative to that of the deepest adsorption site As 
put to 0. Two main low potential trenches are evident along the [110] direction, at both 
sides of the As in-dimer row. The deepest adsorption site Ag is located in the trench near 
the In vacancy of the surface layer and is bonded to the two As atoms of the in-dimer and 



to one In atom (see Fig. 4(a)). Other deep minima (A4, Ag, and Ag) are located along 
the two deepest trenches and, together with Ag, represent the lowest energy set of minima. 
Their adsorption configuration is similar to that of Ag. A second class of adsorption sites 
(Ai, A3, A5, All) is characterized by the In adatom bonded to two other In atoms (see the 
configuration at minimum in Fig. |4(b)D . At the minima Ay and A2 the adatom is bonded 



to four In atoms, reproducing the bonding scheme typical of bulk Indium (Fig. 4(c)). In 
these last cases the corresponding surface electronic structure is metallic. To the last class 
of adsorption sites belongs the shallow minimum Aio where the adatom is bonded to the 



two As of the ad-dimer (see Fig. 4(d)). This configuration plays the role of a precursor for 
the much stronger bonding of In into the ad-dimer (see below). 

We have accurately checked that the shallow minima in the PES are actually minima and 
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FIG. 3: (Color online) (a) PES of an In adatom on the 02 (2 x 4) surface. The PES is drawn for 
4 surface unit cells and a top view of the atomic structure is also shown. The minima are labeled 
Ai. (b) Side view of the surface layer. 

not saddle points of the energy landscape. 



The location of the calculated barriers T„ is shown in Fig. 5(a) The barriers are lower 



along the [110] direction (the highest barrier being 430 meV and the lowest being 52 meV) 
than along the orthogonal [110] direction (see Table [ll|) . This difference is mainly ascribed to 
the presence of the As ad-dimer oriented along the [110] direction, that creates an unstable 
region for the In adsorption, as can be seen in the PES of Fig |3} From these considerations 
we can preliminarily conclude that adatom surface diffusion should be strongly anisotropic 
with the [llO] direction along the trenches favored. 



2. /32(2 X 4) surface reconstruction 

Since the symmetric part of the (32 has the same structure of half of the 02 surface, the 
coordinates of the adsorption sites can be roughly estimated, and their exact location can 
be found with a further relaxation of the surface-adsorbate system. 

Also for the (32 surface reconstruction, we have found 11 adsorption sites for the Indium 



adatom (see Table III and Fig. 5(b)). However, owing to the symmetry properties of the 
system, those minima reduce to only 6 inequivalent minima. 

We first notice that the adsorption sites Ai, A2, A4, Ag, A7, Ag, Ag, are in strict corre- 
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(a) 



(b) 




(c) (d) 

FIG. 4: (Color online) Adsorption configurations for the In adatom in different adsorption sites on 
the a2 surface reconstruction, (a) Configuration (b) Configuration (c) Configuration A2 (d) 
Configuration Aiq 



spondence to those of the ^2 surface reconstruction (see Fig [5]), in the region where the two 
structures are similar i.e. far from the (32 second ad-dimer. The positions of the minima are 
only slightly affected by the presence of the extra As ad-dimer, and the deepest adsorption 
site Ag lies in the same position as for the 02 reconstruction. The minima A3 and A5 are 
shifted towards the trench with respect to the corresponding ones of the 02 reconstruction. 
This is due to the fact that the In adatoms in A3 and A5 on the 02 reconstruction are bonded 
to the buckled In at the trench edge. The same In is not buckled in f32 due to its bond to 



the second As ad-dimer row (compare Fig. 4(b) to Fig. 6(a)). To quantify the difference 



between the A1-A9 adsorption sites in the two reconstructions, we give in Table [IV] the ad- 
sorption energy difference AE = E^{Ai) — E°'{Ai) (negative if In is more strongly bonded 
in the P2 than in 02 reconstruction) and the displacement Ar = yj^ii^i ~ ^f) between the 



corresponding binding locations. The adsorption energy is defined as 



E"(A,) = E{a2 + In^J - E{a2) - ^(In) . 



(2) 



In /32, A3/ FiglSlis the symmetric of the minimum A3. This minimum is not present in the 
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FIG. 5: (Color online) Schematic view of the adsorption sites locations (only i is written) on the: 
(a) 02 (2 X 4) and (b) /?2 (2 x 4) reconstructions. The saddle points are labeled as Tjj. 

0^2 reconstruction because of the lower symmetry. A new minimum A12 is now present for j32 
on top of the two ad-dimers and it is actually quite stable, since the In-adatom is fourfold 



bonded to the surface (see fig. 6(b)). In this case the confinement barriers are calculated to 
be of the order of 135 meV. The configuration for In in the adsorption sites different from 
those of the ^2 reconstruction, are given in Fig. [6j 

Also in the /?2 case, two trenches in the PES are evident along the (110) direction at 
both sides of the As in-dimer. The existence of the A12 adsorption site on /?2 could originate 
additional diffusion paths along the same direction connecting the minimum Au to the two 
sites A2 and A-j. The corresponding saddle point energies are reported in Table |Vj We can 
easily see that the barriers along [110], are in the range 100-580 meV and are lower than 
the barriers along the [110] direction that span the range 100-760 meV. This consideration 
strongly indicates that also for the j32 surface reconstruction, the In-adatom diffusion should 
be highly favored along the (110) direction. 



C. Additional adsorption sites due to In insertion in the As dimers 

We know from the literature [121120] that additional minima could be present, correspond- 
ing to the In adatom breaking the bond between the As adatoms in a dimer and binding 
itself to them. This process requires overcoming an energy barrier, thus its rate should 
depend on the growth temperature. These sites can be of significant interest, since they 
introduce a local modification of the PES that can be very important for surface diffusion. 



10 




(a) 



(b) 



FIG. 6: (Color online) Adsorptions configurations for the In adatom on the /?2 surface reconstruc- 
tion, (a) Configuration (b) Configuration An 

For example, such sites could be in a bridge location between two other adsorption sites, 
thus introducing a further minimum in an otherwise thicker barrier. Such an additional 
adsorption site can affect the diffusion kinetics of the In adatom on the surface. 

We have considered both the As in-dimers and ad-dimers of the 02 and (32 surface re- 
constructions and determined the energy necessary to break the As-As bond and bind the 
In adatom. First of all we have put In in the middle between the two As atoms, in a 
vertical plane containing all the three atoms. We have then relaxed the system obtaining 
the minimum energy position for the In adatom. Finally we have performed different NEB 
calculations in order to find the energy barriers between these adsorption sites and the neigh- 
boring sites. We have found that the adsorption configuration with the In atom bonded to 
the two arsenic of the dimers are in general more stable than most of the other adsorption 
sites (see Tables [l] and III). 



1. 02(2 X 4) surface reconstruction 

In this case we have found 2 additional stable sites for the In adsorbate. The first one 
corresponds to the breaking of the As ad-dimer and the second one to the breaking of the 
As in-dimer. The two cases are shown in figures [T] and |8| respectively. Figures [7]^a) and 
[8|^a) show the location of these minima, indicated as A^ and Aj respectively, on the surface 
(a =ad-dimer i =in-dimer). Figures [Tj^b) and |8](b) show the atomic configuration of the 
adsorbate bonded to the As atoms in the dimers. The transitions to/from neighboring 
adsorption sites are also drawn in the figures (Figs, ^c) and|8](c)) and the corresponding 
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FIG. 7: (Color online) Additional adsorption site Aa due to the In insertion into the As ad-dimer 
on the 02 reconstruction, (a) The location of the adsorption site is indicated by a star, (b) Atomic 
configuration, (c) Scheme of the transitions between the A^ adsorption site and the neighboring 
sites, (d) Energy profile for barrier Tas from the precursor site Aio to Aq. 

minima and saddle point energies are given in Tables [T] and [TT] respectively. We can see 
that those minima are deep and stable, since the confining barriers (~ 0.8 eV), pictured 
in Figures and [Sj^d) are high. The energy barrier needed to fall into from the site 
Aio (the precursor site) is not too high (around 200 meV) and this process is likely to occur 
even at relatively low temperature. 



2. /32(2 X 4) surface reconstruction 

In this case we have two As ad-dimers instead of one. We have found one adsorption site 
for the In adatom for each one of the ad-dimers indicated in Figures [9] as A^i and Aa2- We 
have found that the In adsorbate does not lie in the vertical plane passing through the As 
dimers, but the bonds are tilted towards the other ad-dimer, that is towards the symmetry 
plane passing through the A12 adsorption site. Regarding the As in-dimer we have found 



two distinct possibilities for In adsorption that we indicate as Aji and Aj2 in Figure 10 As 



can be seen from Fig. 10 b), the In atom bonds to the two As of the in-dimer forming a 
tilted angle with respect to the vertical z direction. Since the (110) is a plane of symmetry 
for the system, the two configurations for the adsorbate, corresponding to angles (p and —(f), 
with respect to the z axis are equivalent and have the same probabiUty to occur. 
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FIG. 8: (Color online) Additional adsorption site Aj due to the In insertion into the As in-dimer 
on the 02 surface reconstruction, (a) The location of the adsorption site is indicated by a star, (c) 
Atomic configuration, (b) Scheme of the transitions between the Aj adsorption site and the neigh- 
boring sites, (d) Energy profile for barriers Tn and Ti2 from the two besides in-dimer adsorption 
sites A4 and Ag to Aj. 



Figures [9[a) and[lO[a) show the location of the adsorption sites for the As ad-dimers and 
As in-dimer, respectively. The atomic configurations for the adsorbates are given in Figures 
[oj^b) and[IO|^b), whereas the transitions to/from neighboring adsorption sites are reported in 
Figures [9|c) and [To|^c). Finally Figures |9](d) and 10 d) plot the energy barriers separating 
these adsorption sites from the neighboring ones. 



Tables III and |V] report the energies for the stable configurations and the corresponding 
saddle points. Even in this case we have found very high binding energies for these adsorption 
sites and high energy barriers (around 0.7 eV) that make these configurations very stable 
once the adatom reaches them overcoming the barriers. 



IV. DISCUSSION 



In order to evaluate the effects of the strain and the role of the different surface recon- 
structions on the In adsorption properties we now compare our results with the results of 
similar calculations reported in the literature. To better evaluate the differences we report in 



Fig. 11 the more common adsorption features as extracted by the literature and the present 
work. In the same figure we name each specific adsorption site with respect to its position 
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FIG. 9: (Color online) Additional adsorption sites A^i and A^i due to the In insertion into the As 
ad-dimer on the /32 surface reconstruction, (a) The location of the adsorption sites are indicated by 
stars, (b) Atomic configuration, (c) Scheme of the transitions between the A^i and A^i adsorption 
site and the neighboring sites, (d) Energy profile for barriers Ta2 from the top ad-dimer adsorption 
site Ai2 to Aa- 



relative to the As in-dimers (i) and ad-dimers (a), and indicate also the reconstructions 
where these adsorption sites are typical. 



A. 02 vs (32 reconstruction (this work) 

There are no big differences between the adsorption features on the a2 and on the /32 
surface reconstructions. In particular, the adsorption energy and location differences (see 



Table IV) between corresponding adsorption sites, are within a few hundreds of meV. The 
greatest differences are those in the region of the missing ad-dimer where the structural 
difference between the two surface reconstructions is larger. On the (32 reconstruction the 
positions and energies of the adsorption sites are modified owing to the symmetry properties 
of the reconstruction. The most important difference is the presence of the adsorption sites 
A3/ and A12 on the (32 surface reconstruction, that are absent on the 02 surface reconstruction. 
On the contrary, the adsorption sites Aio and An of 02 disappear in P2 transforming into 
the top ad-dimer site A12. 
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FIG. 10: (Color online) Additional adsorption sites Aji and Aj2 due to In insertion into the As 
in-dimer on the P2 surface reconstruction, (a) The locations of the adsorption sites are indicated 
by stars, (b) Atomic configuration, (c) Scheme of the transitions between the Aji and Aj2 adsorp- 
tion sites and the neighboring sites, (d) Energy profile for barrier Tn from the besides in-dimer 
adsorption site A4 to Ai2. 

B. 1.75 ML InAs on GaAs (001)/32(2 x 4) vs InAs (001)/32(2 x 4) 

In this subsection we compare our results for In adsorption on the /32(2 x 4) reconstruction 
of the InAs WL on GaAs(OOl) with the results for In adsorption on the (32{2 x 4) recon- 
struction of pure InAs(OOl) of ref. [211 E2|- The main difference between the two cases is 
that the InAs WL on GaAs is subjected to a compressive strain on the (001) plane. Thus, 
the comparison allows us to evaluate the effect of the WL strain on In adsorption. We 
principally refer to Fig.2 of ref.|2I] and to Table VI of the present paper, where we have 



reported the energies of the minima of In on the InAs surface 

In the case of pure InAspI], only six minima were calculated and their positions on the 
surface are almost the same as those calculated by us. In particular, only one minimum 
Bg in the bridge in-dimer position (see Fig. [II] ) was found halfway between the two equal 
minima Ag and Ag calculated by us for the InAs/GaAs case. The two minima Ag and Ag 
are separated by a relatively high barrier of 300 meV in our calculation. We recall here that 
for a2 we calculate barriers only 42 and 142 meV high between the same adsorption sites. 
(Ag and Ag are no more equivalent in 0:2). 

In the case of pure InAs, the besides in-dimer B4 and Eg show the strongest binding. In 
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FIG. 11: Scheme of the main adsorption sites for an In adatom on (001) reconstructed III-V 
surfaces: (a) besides in-dimer and bridge in-dimer locations, (b) top ad-dimer location, (c) bridge 
ad-dimer location, and (d) missing ad-dimer locations. 



the case of the InAs WL on GaAs, instead, the more strongly bonded sites are located at 



the bridge in-dimer positions Aq and Ag. From Tables |III| and |VI| we see that the energy 
difference between the sites (4-9) and (6-8) is larger for the InAs WL on GaAs case (240 
meV vs 61 meV). It is reasonable to think that this difference is to be ascribed to the strain 
originated by the lattice mismatch between InAs and GaAs. Also the absence of the shallow 
adsorption sites corresponding to Ai, A5, A3 and Aio could be originated by the absence of 
the strain. 

Nevertheless, it is possible that these differences are originated also by the different grid 
resolution used for calculating the PES, or by the interpolation algorithms that sometimes 
are not able to identify correctly the minima or by the difference between GGA used in ref. 
[2T] and LDA used in this work. 
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C. In on InAs/GaAs /?2(2 x 4) vs Ga on GaAs /32(2 x 4) 



At this stage of the discussion we compare our results for the j32 reconstruction with 
the resuhs obtained in ref. [2U], where the PES of a Ga adatom on the /52(2 x 4) surface 
reconstruction of GaAs(OOl) was studied. This comparison evidences the different behavior 
of In and Ga adsorbates, as opposed to the WL composition. 

In the case of Ga on GaAs (001) only three minima were reported (without considering 
the additional minima where Ga breaks the dimers). Two minima are at bridge ad-dimer 
positions while the third one is at a bridge in-dimer position. These adsorption sites corre- 
spond to the A2, A7 and Ag-Ag minima on the InAs WL and to B2, By and Eg minima of 
the In on InAs(OOl) case |2T]. These minima are quite stable. Notice the missing of the A12 
adsorption site on the top ad-dimer location, which is present in the case of the In adsorbate 
on /?2 reconstructed InAs covered surfaces. Also the adsorption sites in the besides in-dimer 
locations are missing in this case. 



D. In adatom on InGaAs surfaces 

In this section we compare with the results of ref. [121 [IB] where InGaAs reconstructions of 
the WL on GaAs were considered. The PESs of the In adatom on: (i) In2/3Gai/3As(001) (2 x 
3) (corresponding to a In coverage 9 = 0.66), (ii) In2/3Gai/3As(001) (1 x 3) (corresponding to 
6 = 0.66), and (iii) GaAs(OOl) c(4 x 4) (corresponding to ^ = 0) surface reconstructions were 
studied. From all these calculations we can see that one adsorption site is always located 
in a top ad-dimer position, in the same configuration of the A12 adsorption site on of 
this work. The depth of these adsorption sites defined by the height of the lowest barrier 
confining the adsorbate are: 140 meV, 60 meV, 100 meV, 385 meV and 135 meV for (2x3), 
(1 X 3), c(4 X 4), InAs /52(2 x 4) and InAs WL on GaAs /?2(2 x 4), respectively. We can see 
that these minima are not strongly bonded, apart from the case of In on pure InAs(OOl) (32 
reconstructed surface. 

The deepest minima are located at the besides in-dimer location, both in the case of the 
(2 X 3) reconstruction of [12] and of the reconstructions studied in this work. In the case 
of the c(4 X 4) reconstruction, stable locations are at the missing ad-dimer sites (see fig. 



11). In these adsorption sites the In adatom binds to four As atoms below. In the case of 
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the (1x3) reconstruction the stable location sees the In atoms bonded to the cations layer 
below. 

V. CONCLUSIONS 

We have carried out a first-principles study of the adsorption properties of an In adatom 
deposited on a InAs wetting layer on GaAs(OOl) reconstructed 02(2 x 4) or /32(2 x 4). We 
have first studied the equilibrium geometry of the surfaces and the relative phase stability 
with respect to the growth conditions. We have found that, in the case of the InAs WL 
the interlayer distances are systematically higher than those of pure GaAs, owing to the 
compressive strain. Moreover, we have found that the a2 reconstruction of the InAs WL is 
favored over the fi^ reconstruction, for a wider range of growth conditions, when compared 
to the bare GaAs(OOl) surface. 

We have analyzed the PES for an In adatom determining the adsorption sites and the 
saddle points for both 0.2 and ^2 reconstructions. We have found 11 minima in the PES 
for both reconstructions. Most of them occur in the same positions in both cases. Two 
minima are instead specific of the a2 and reconstructions, due to the presence on of 
an additional As ad-dimer. For the adsorption sites present on both reconstructions, we 
have studied the change of the adsorption energy and location caused by the presence or 
absence of the second As ad-dimer. We have identified two low-energy trenches along the 
[110] direction, alongside the in-dimer chain, for both 0:2 and /32- 

The comparison to analogous studies reported in the literature has revealed some general 
features for the In adsorption on different surface reconstructions. In particular, we have 
found that stable adsorption sites arc always located: (i) besides the in-dimers, (ii) at the 
bridge positions between in-dimers, (iii) between ad-dimers (only for ^^2 which has two ad- 
jacent As ad-dimers), and (iv) at the bridge between ad-dimer positions. We have identified 
also other shallower adsorption sites which are more reconstruction specific and are related 
to a2 having a lower symmetry than ^2- 
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site energy (meV) 



Ai 


220 


A2 


387 


A3 


312 


A4 


49 


A5 


344 


Ae 


99 


A7 


559 


As 





Ag 


112 


Aio 


868 


All 


498 


A. 


-100 




239 



TABLE I: Adsorption energies for the In adatom in the different adsorption sites of the 02 sur- 
face reconstruction, relative to the As energy, labeled as A^. The labels A^ and A^j refer to the 
configurations of the adsorbate into the in-dimer and the ad-dimer, respectively. 
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APPENDIX A: ENERGIES OF MINIMA AND SADDLE POINTS 

In this appendix we report all the numerical values calculated for minima and saddle 
point energies. 
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saddle 


energy (meV) 


T2,9 


895 


T2J 


676 


T3,7 


635 


T3,4 


456 


T4 9 


498 


^4,5 


453 


Ti 10 


893 


Tg^g 


142 


T8,9 


428 


Tl,9 


272 


Tl,8 


418 


T5 


554 


Tio.ii 


929 


Ty 11 


743 


Til, 7 


750 


Til 


613 


Ti2 


493 


Tal 


893 


Ta2 


926 


Ta3 


1117 



TABLE II: Energies of the saddle points in the PES of the 02 surface reconstruction, relative to 
the As adsorption energy. The saddle points connecting the minima An and Am are labeled by 
Tn,m) while the labels Tj„ (Tan) indicate the saddle points between the minima into the in-dimer 
(ad-dimer) and the neighboring sites, as shown in Fig. [s] (Fig. [7| 
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site 


energy (meV) 


Ai, As 


216 


A2, Ay 


506 


M, A-j/ 

J 


323 


A4, Ag 


240 


Ae, Ag 





A12 


588 




-330 


A„ 


158 



TABLE III: Adsorption energies for tlie In adatom in the different adsorption sites of the P2 
surface reconstruction, relative to the As energy, labeled as A„. The labels Aj and Aq refer to the 
configuration of the adsorbate into the in-dimer and the ad-dimer, respectively. 



site 


AE (meV) 


Ar/ao 


Ai 


-134 


0.022 


A2 


-7 


0.069 


A3 


-121 


0.330 


A4 


66 


0.168 


A5 


-259 


0.350 


Ae 


-224 


0.100 


A7 


-179 


0.588 


As 


-125 


0.018 


A9 


4 


0.031 



TABLE IV: Adsorption energy differences and position displacements between the In adsorbate 
for the corresponding adsorption sites of the P2 and the 02 reconstructions. 
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saddle 


energy (meV) 


T4,5 


392 


^5,6 


424 


Tl,l2 


976 


Tl2,2 


723 


T2 11 


727 


^2,3' 


974 


^2,7 


589 


Tg 3 


577 


^4,3 


424 


T4 9 


701 


T6,8 


296 


Til 


589 


Ti2 


94 


Tal 


760 


Ta2 


610 


Ta3 


1024 



TABLE V: Energies of the saddle points in the PES of the (32 surface reconstruction, relative to the 
As {(32) adsorption energy. The saddle points connecting the minima An and Am are labeled by 
Tn,m) while the labels Tj„ (Tan) refer to the saddle points between the minima into the in-dimer 



(ad-dimer) and the neighboring sites, as shown in Fig. 10 (Fig. 



site energy (meV) 
B2, B7 302 

B4, Bg 

Bg 61 
Bii 163 

TABLE VL Adsorption energies (calculated with respect to B4 energy) for the In adatom in the 
different adsorption sites on the InAs (32 surface reconstruction, labeled as B„ (see ref. |2T] [33]). 
The positions of the minima B„ corresponds to the minima A„ in fig. |5(b) 
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